ENUMERATION SCHEMES AND, MORE IMPORTANTLY, THEIR AUTOMATIC GENERATION 



Doron ZEILBERGER^ 

Towards the Virtual Combinatorialist 

'oO ! is soon to teach our computers how to become full-fledged humans. It is even premature 

0^ ' to teach them how to become mathematicians, it is even unwise, at present, to teach them how to 

become combinatorialists. But the time is ripe to teach them how to become experts in a suitably 
defined and narrowly focused subarea of combinatorics. In this article, I will describe my efforts to 
^ ■ teach my beloved computer, Shalosh B. Ekhad, how to be an enumerator of Will classes. 

' Just like the Wright Brothers' first flight, and Len Adelman's first DNA computer (and perhaps like 

■ the first Fleischmann-Pons Cold Fusion experiment), the actual accomplishments are very meager. 

But I do hope to demonstrate feasibility, and jump-start a research area that is not quite part 
of AI, and not quite algorithmic combinatorics. It is something brand-new, let's call it Artificial 
Combinatorics (AC). The closest analogy is Deep Blue, but instead of playing Chess, Shalosh will 
play a game called 'Wilf-class enumeration'. It is still very far from beating the Kasparovs of the 
area (Miklos Bona, Alex Burstein, John Noonan, Prank Schmidt, Rodica Simion, Zved Stankova, 
Julian West, to name a few), but I am sure that the ultimate version will. 
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With all due respect to Wilf classes and Enumeration, and even to Combinatorics, the main point 



^ . of this article is not to enhance our understanding of Wilf classes, but to illustrate how much (if 

i not all) of mathematical research will be conducted in a few years. It goes as follow. Have a (as 

•/^ I of now, human) mathematician get a brilliant idea. Teach that idea to a computer, and let the 

^ I computer 'do research' using that idea. 

0^ ' 

fH , Accompanying Software 



In a sense, this article is a user's manual for the Maple package WILF, available from my homepage 



http : //www .math. temple . edu/~ zeilberg (click on Maple packages and programs, and then 



download WILF). The empirical version of the rigorous WILF is HERB, also available there. 



S 
> 

X 

• Enumeration Schemes 

Suppose that we have to find a 'formula' (in the sense of Wilf[Wi], i.e. a polynomial-time algorithm) 
for computing a„ := where An is an infinite family of finite sets, parameterized by n. Usually 
An is a natural subset of a larger set and is defined as the set of members of Bn that satisfy 
a certain set of conditions C„. For example if An is the set of permutations on {1,2, . . . ,n}, then 
Bn may be taken as the set of words of length n over the alphabet {1, 2, . . . , n}, and C„ can be 
taken as the condition: 'no letter can appear twice'. A naive algorithm for enumerating An would 
be to actually construct the set, by examining the members of Bn, one by one, checking whether 
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they satisfy C„, and admitting those that quahfy. Then a„=Cardinality of A^. 

But a much better approach would be to find a structure theorem that expresses An, using unions, 
Cartesian products, and possibly complements, of well known sets. Failing this, it would be also nice 
to express A^, recursively, in terms of A^-i, An-2, ■ ■ ■, and easy-to-count sets, getting a recurrence 
formula. Going back to the permutation example, Levi Ben Gerson ([L]) proved the structure 
theorem An = {1, 2, . . . , n} x An-i, from which he deduced the recurrence an = na„_i, enabling a 
polynomial- (in fact linear-) time algorithm for computing Oj, for 1 < i < n. 

Alas, this is not always easy, and for many enumeration sequences, e.g. the number of self-avoiding 
walks, may well be impossible, and who knows, perhaps one day even provably impossible. 

It is conceivable, however, that a combinatorial family A(n), does not possess a recursive structure 
by itself, but by refining it, using a suitable parameter, one can partition A{n) into the disjoint 
union: 

n 

Ain) = \jB{n,i) , 
1=1 

and try and find a structure theorem for the two-parameter family B{n,i). This will imply a 
recurrence for the cardinalities h{n,i) := \B{n,i)\, that would enable a fast algorithm for a(n) = 

Sometimes, not even the B(n,i) suffice. Then we could try to partition B{n,i) into the following 
disjoint union: 

i — \ n 

B{n,i) = \JCi{n,i,j)U (j C2{n,i,j) , 

j=l j=i+i 

and try to express Ci(n,i,j) andC2(n,z, j) in terms of ^(m), B{m,i'), Ci{m,i' and C2(m, 
with m < n. One can keep going indefinitely. If this process halts after a finite number of refine- 
ments, then we have indeed a formula (in the sense of Wilf) for a{n). 

Permutations with Forbidden Patterns; Wilf classes 

A pattern of length A; is a permutation of {1,2,..., k}. 

The reduction of a vector of distinct integers of length k is the pattern obtained by replacing the 
smallest element by 1, the second smallest by 2, . . ., the largest by k. For example, the reduction 
of 264 is 132. 

A permutation vr of length n, of an ordered set, is said to contain the pattern a of length k if there 
are places ii < . . . < ik, such that ni^TTi^ . . . TTj^. reduces to a. For example 51872463 contains the 
pattern 3421 (with ii = l,i2 = 3, i^ = 6,ii = 8). 

Our Goal is to investigate the following 

Problem: Given a set of patterns V, study the sequence a{n;V) := |^(n;7^)|, where A{n;V) is 
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the set of permutations on {1,2,..., n} that avoid the patterns of V. 



We will sometimes write A(n) and a(n) instead of A{n;V) and a{n;V), where the V is implied 
from the context. 

Of course a(n;0) = n\, and a(n; {12}) = 1. It is well known (see [We][K]) that a(n; {123}) = 
a(n; {132}) = C„, where C„ = (2n)!/(n!(n + 1)! are the Catalan numbers. 

Naively, it takes an exponential time to compute the first n terms of the sequence a{n;V). The 
best possible scenario is an explicit formula for the sequence. Failing this, it would be almost as 
nice to have a linear recurrence equation with polynomial coefficients (and hence inducting it in 
the hall of fame of P-recursive sequences) . 

Sometimes it happens that there exist two sets of patterns V and V' such that a{n;V) = a{n;P') 
for all n. If this happens for a non-trivial reason, we say that V and V belong to the same Wilf 
class. 

The approach that I am going to investigate here, and that I taught Shalosh B. Ekhad, is to 
partition the set A{n; V) into a finite disjoint union of subsets, indexed by certain prefixes, and to 
try to deduce rigorously, certain recurrence relations between them, that would enable one to set 
up an enumeration scheme, that would lead to a polynomial-time algorithm for computing a{n;'P). 
In the future it is hoped that these schemes could also be used to derive (once again rigorously) 
generating functions for the more refined quantities, that should lead to a constant-term expression 
for a(n; V), and that would make it decidable whether V and V belong to the same Wilf class. 

Apology: The success rate of the present method, in its present state, is somewhat disappointing. 
Ekhad was able to reproduce the classical cases, and a few new ones, but for most patterns and 
sets of patterns, it failed to find a scheme (defined below) of a reasonable depth. But the present 
set-up for setting up a scheme could be modified and extended in various ways. We do believe that 
an appropriate enhancement of the present method would yield, if not a one-hundred-per-cents 
success rate, at least close to it. 

The Most Trivial Non- Trivial Example 

In order to illustrate the method, let's work out in full detail, and in plain English, an Enumeration 
Scheme for the set of permutations avoiding the pattern 123. 

Let A{n) be the set of such permutations. Break A{n) up into the union 

A(n) =ur=iAi(n;z) , 
where Ai{n,i) are the set of 123-avoiding permutations on [l,n] that start with i, i.e. 

Ai(n,?) := {tt G A(n);7ri = i} . 
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Let's break up Ai{n;i) into the following disjoint union: 



i-l 



Ai{n,i) = [J A2i{n;h,i)U [J Ai2{n;i,j) , 

h=l j=i+l 

where for 1 < i < j < n, 

Ai2{n,i,j) := {tt € A{n);ni =1,7x2= j} , 
A2i{n,i,j) := {tt G A(n);7ri = j , 172 = i} . 

Now let's look at a typical element tt of A2i{n,i, j). It has vri = j and 7r2 = i. Deleting the first 
entry (i.e. the j), and reducing, gives us a member a of Ai{n — 1, i), since you can't create trouble 
(i.e. a 123 pattern) by skipping town. FURTHERMORE, if you start out with any element a 
oi Ai{n — make room for the j (by raising all members > j by 1), and insert j at the very 
beginning, you can't create new trouble, because the j is "body-guarded" by the i. Let me explain. 
If inserting j at the very beginning would have caused an undesirable 123-pattern, then of course, 
the "1" of the "123" would be j, but the "2" can't be i (since j > i), hence in this hypothetical 
new 123 that the insertion of j at the very beginning would have created, the "2" would have been 
a certain aa and the "3" would be a certain at,, with 2 < a < b. But had that been the case, then 
i would have been an EVEN BETTER "1"!, because i < j, and ai < aa < ai, would have already 
been a 123-pattern. It follows that A2i{n;i,j) is isomorphic to Ai{n — 

Now let's examine Ai2{n;i, j). If j < n, then n must be somewhere to the right, making ijn a 
123-pattern. Hence Ai2{n;i,j) is empty when j < n. On the other hand, if j = n, then deleting it 
creates a member oi Ai{n — Conversely, given any member of Ai(n — l;i), inserting n right 
after the first entry (hence making n the second entry), can't create a 123. Hence Ai2{'n;i,n) = 
Ai{n — Summarizing, we have the following structure: 

n 

A{n) = [jAi{n;i) , 

i=l 

i — 1 n 

Ai{n,i) := [J A2i{n;h,i) U [J Ai2{n;i,j) , 

h=l i=i+l 

A2i{n;i,j) = Ai{n-l;i) , 



Ai2{n;i,j) = 



0, if j < n; 

Ai{n-l;i), ifj = n. 



Now let a(n) := \A{n)\, ai(n;i) := \Ai{n;i)\, ai2{n;i,j) := \Ai2{n;i,j)\, a2i{n;i,j) := \A2i{n;i,j)\. 
Taking cardinalities gives us the following scheme: 

n 
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i— 1 n 

ai{n,i) ■.= '^a2i{n;h,i) + ^ ai2(n;i,j) , 

h=l j=i+l 



02i(n;i,j) = ai(n - 
ai2(ra;i,j) = 

Since 021 and ai2 can be expressed in terms of ai , the above can be condensed to 



0, if j < n; 

ai(n — 1; i), if j = n. 



i=l 

\ i 

ai(n, i) = I ai(n — 1; /i) j + ai(n — 1; i) = ai(n — 1 ; /i) , 

\h=l J h=l 

from which it follows that ai(n;i) — ai{n;i — 1) = ai(n — l;z). This recurrence and the obvious 
boundary conditions, ai(n;l) = 1 and ai(n;n+ 1) = 0, yield the explicit solution ai(n;i) = 
- ("+r') Hence a{n) = a,{n + l,n + 1) = = (^:)/(n + 1). 

Introducing Prefix Schemes 

The above example leads to the following definitions. Fix a set of patterns V. Let o" = o"i . . . cr^ be 
any permutation of length k. For 1 < ii < ^2 < . . . < ifc < n, let 

A„{n;V;ii, . . . ,ik) ■= {tt G A{n;V) \ tti = Vi,7r2 = ia^, ■■■ , T^k = • 
For example ^132(6; {1234, 1432}; 2, 3, 5) := {25314, 25341}. 

The set of refinements of a permutation a = ui . . .Uk of lengths A; is the set of A; + 1 permu- 
tations of length k + 1 that have the property that deleting their last entry reduces to a. We 
will denote the set of refinements of a by Refinements(a). For example, Refinements{312) = 
{3124,4123,4132,4231}. 

For any set of patterns V, and any prefix-permutation a, we have the obvious recurrence 

fe+i V— 1 

A^{n;V;ii,...,ik) = [j [j A^^r-){n;V;ii, . . . ,ir-i,j,ir, ■ ■ ■ ,ik) , (*) 

r=l j=ir-i+l 

where crW is that element of Refinements{a) that ends with r. We agree that iq = and ik+i = n. 
For example 

n — 1 12—1 *3 — 1 

Ai32{n-ii,i2,i3)=[jA243i{'n;j,ii,i2,i3)^ [J Aii32{n;ii, j, 12,^3) ^ |J ^i423(?i;n,«2, j,«3) U 

n 

[J Ai324:{n;V;ii,i2,i3,j) ■ 

j=i3+l 



Given a set of patterns, V, and a prefix permutation cr, it may happen that all the members of 
Aa{n; V), for all n, must have either one or more of the properties ii = 1, or ii = 12, 12 = h, ■ ■ ■ ,ik = 
n. Let J(cr) be the set of < r < /c such that v = V+i is forced for all members of A^{n; V). In the 
above example with V = {123}, we had J(12) = {2}, since every 123-avoiding permutation that 
starts with ^1^2, where ii < 12, must have ^2 = n. Another Example: If P = {1234, 1324, 1243}, 
then J(2413) = {4}, since every member of A24i3{n;T';ii,i2,i3,i4) must have its second entry, 14 
be equal to n, or else the n would be to the right of the fourth entry, that would cause 12141311, to 
be an illegal 1324. 

Consider all permutations. Let's remember our goal. Given a set of patterns V, we want to count 
the number of permutations that avoid all the patterns of V. Let's call these permutations law- 
abiding. Hence a permutation that has one or more occurrences of the patterns V should be called 
criminal. Obviously, deleting any entry of a law-abiding permutation (and reducing) , gives rise to 
another law-abiding permutation. Hence given a prefix permutation a = ai . . .ak, then for each 
place r (1 < r < fc), we have the obvious inclusion: 

Fr{A^{n;V;ii,. . . ,ik)) C A^^^^{n - l;V;ii, . . . , v+i - 1, . . . ,Zfc - 1) , {Rodica) 

where Fj. is the injective mapping that consists of deleting (wherever it is), and reducing, and 
(T(t.) is the permutation obtained from a by deleting r and reducing. 

BUT, it may happen that the inclusion {Rodica) is an equality. This leads to the following 

IMPORTANT DEFINITION. Given a set of patterns V, and a prefix permutation a, the place 

r (1 < r < fc) is reversely deleteable if for all 1 < ii < . . . < < n, 

Fr{Aa(n\V\ii,. . .,ik)) = Ac,^^^{n - l;V;ii,. . . , V-i,ir+i - 1, • • • ,^fe - 1) • (Julian) 

We must find a way to know, a priori, by using logical reasoning, whether a place is reversely 
deleteable. 

How can we be sure that inserting v in any permutation of A(j^^^{n — 1; "P; ii , . . . , v+i — 
1, . . . , ife — 1), in the appropriate place, is a safe thing to do? We must be assured that its insertion 
cannot cause any trouble. How can we be sure? We have to look at all the conceivable events 
that its insertion can cause, i.e. all the possible forbidden patterns of V, with one of its entries 
coinciding with the inserted v- If for each such possible scenario, in which v participates, we can 
logically deduce the existence of another event, in which v does not participate, then we are safe, 
since we are covered. 

Example: Let V = {123}. For a = 21, the place r = 1 is reversely deleteable, since suppose that 
inserting the 12 would have created a 123. Then there exist a and 6 in tt such that i2ab is a 123- 
pattern. Had that been the case, then kal vakhomer, i-^ab is also a 123-pattern. 

Another Example: Let V = {1234, 1324, 1243}. a = 2413. We saw above that J(c7) = {4}, i.e. 
the permutations tt that we are examining start with i2niii3 for some 1 < ii < i2 < 13 < n. 1 claim 
that the first place, i.e. the Z2 is reversely deleteable. 
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Let's look at the kind of trouble the insertion of 12 at the beginning could cause. 

Event 1: There exist 4<6<c<d<ra such that ^27^a7^b7^c is a delinquent 1234. Then 12 can be 
bailed out by ii, since iiTTairbT^c would have been a bad event even before the insertion. 

Event 2: There exist 4 < 6 < c < n such that i2i?,T^b'^c is a delinquent 1234. Then i2 can be bailed 
out by ii, since iiisTTbTTc would have been a bad event even before the insertion. 

Event 3: There exist A<a<b<c<n such that ^27^a7^67^c is a delinquent 1324. Then 12 can be 
bailed out by ii, since ziTTiTTftTTc would have been a bad event even before the insertion. 

Event 4: There exist 4 < 6 < c < n such that i2hT^b'^c is a delinquent 1324. Then Z2 can be bailed 
out by zi, since iiizi^bT^c would have been a bad event even before the insertion. 

Event 5: There exist 4<a<6<c<n such that i2'^a'^b'^c is a delinquent 1243. Then {2 can be 
bailed out by ii, since iiTTi-KbT^c would have been a bad event even before the insertion. 

Event 6: There exist A < b < c < n such that i2«37i"67i"c is a delinquent 1243. Then 12 can be bailed 
out by ii, since i\izT^bT^c would have been a bad event even before the insertion. 

So any law-abiding permutation tt, that starts with a triplet that reduces to 312 has an a priori 
guarantee that sticking an entry at the very front, that make it have prefix that reduces to 2413, will 
not get the permutation tt in trouble. Hence we can know for sure that, with V = {1234, 1324, 1243}, 

X / 0, if i4 < n; 

A3i2(n - 1 ; P;zi,Z3 - 1,Z4 - 1), \ii^ = n. 

Two 'Trivial' Examples 

Let's rederive Levi Ben Gerson's famous recurrence for a(n;0), i.e. a(n;0) = na{n — 1;0). 
First, let's refine A{n\ 0) into 

^(n;0) = U-=iAi(n;0;zi) . 

Now, note that r = 1 in cr = 1 is reversely deleteable, since inserting an i\ at the front can't cause 
any trouble, since nothing is forbidden. Hence 

Ai(n;0;ii) = ^(n-l;0) , 

from which follows that 

n n 

a{n; 0) = ^ ai(n; 0; ii) = ^ a(n - 1; 0) = no(n - 1; 0) . 

ii=l ii=l 

Next, let's find a scheme when V = {12}, that is let's try to compute a(n; {12}), the number of 
permutations on n objects with no occurrence of the pattern 12, i.e. the number of decreasing 
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permutations. First, as always 



A(n;{12})= U Ai(n;{12};zi) . 
ii=i 

Now J(l) = {2}, since any 12-avoiding permutation must have ii = n, or else iin would be a 
forbidden pattern. Also r = 1 is a reversely deleteable place, since inserting ii = n at the very 
beginning of a permutation in A{n — 1;{12}) can't possibly cause trouble. We hence have the 
following recurrence: 

A(n;{12}) = A(n-l;{12}) , 

which implies that a{n; {12}) = a(n — 1; {12}), and since a(0; {12}) = 1, we have proven, rigorously, 
that a(n; {12}) = 1, for ah n > 0. 

A prefix scheme can be defined independently of sets of forbidden patterns V as follows. 
The Formal Definition of Prefix Scheme 

Definition: A PREFIX SCHEME is a five-tuple [Redu,Expa,A,B,C] such that: 

(i) Redu is a finite set of permutations (of various lengths) . 

(ii) Expa is another finite set of permutations (of various lengths). It must include the empty 
permutation 0. Expa and Redu are disjoint. 

(iii) ^ is a table assigning to each member cr = cri ... crfc of Redu, a place A[a] in a, i.e. an integer 
i between 1 and k. 

(iv) S is a table assigning to each permutation a = ai . . . ak of Expa, its fc + 1 refinements. 

(v) C is a table that assigns to each member a = ai . . . oi ReduU Expa a certain set J{cr), which 
is a (possibly empty) subset of {0, 1,2,..., k}. 

(vi) For each a G Expa, each member B[a] must belong to Expa U Redu. 

For example, the following [Redu, Expa, A, B,C] is a prefix scheme: Redu = {12,21}, Expa = 
{</>, 1}, ^[12] = 1, ^[21] = 1, B[<P] = {1},B[1] = {12, 21}, C[12] = {2}, C[21] = C[l] = C[</>] = 0. 

We are now ready to interface the abstract notion of prefix scheme to that of forbidden patterns. 

Another Important Definition: 

A Prefix Scheme [Redu, Expa, A, B,C] is an Enumeration Scheme for a set of forbidden 
patterns V, if the following conditions hold: 

(i) For every a G Redu, let a' be the reduction of the permutation obtained from a by deleting the 
entry at the place specified by A[a]. The following property holds: If you take any permutation tt 
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that avoids the patterns in V, and that has a prefix that reduces to cr', and insert a new entry right 
before the {A[a] + 1)*'* place of vr, in such a way that the new permutation has a prefix that reduces 
to 0", then that new permutation also avoids the patterns of V. In other words, the insertion is 
always safe. 

(ii) Suppose that a permutation tt of {1, 2, . . . , n}, avoids the patterns of P, and has a prefix that 
reduces to one of the permutations a = ai . . . ak of Expa U Redu, and the first k entries of vr are 
ia-iicr2 ■■■icrk^ whcrc, of course, 1 < ii < i2 < ■ ■ ■ < ik < n. Then if G C[a], we MUST have 
ii = 1, and if A; G C[a], we MUST have ik = n, and for any other member j of C[a], we MUST 
have ij = ij+i. 

Implementing the Scheme 

Once the computer (or, in simple cases, the human) has found a prefix-scheme for a set of patterns 
V, then we have a "formula" (in the Wilfian sense discussed above) for enumerating a{n;V). 

If (T = cji . . . (jfc G Expa, let a^^'' [1 < j < k + 1), he that member of B[a] that ends with j. We 
have, for 1 < ii < . . . < ifc < n 

fc+i v-i 

aa{n;T';ii,...,ik) = ^ ^ a^(j){n;V;ii, . . . ,ij-i,r,ij, . . . ,ik) , 

j=l r=ij-i+l 

where iq = and ik+i = n + 1. 

If cr = (Ti ... (Tfe is in Redu, let A[a] = r. If C[a\ = 0, then 

aa{n;V;h,. ■ ■ ,ik) = aa^r){'^ - • • • , V-i, V+i - 1, ■ ■ ■ ,«fe - 1) , (Miklos) 

where a"(r) is the permutation obtained by reducing the permutation obtained from a by deleting 

If J := C[(t] 7^ then {Miklos) holds whenever . . ,ik) "obeys J", i.e. whenever ij = ij^i for 
all j & J (recall that io = 1, ik+i = n). Otherwise aa{n;V;ii, . . . ,ik) = 0- 

Since the sets Redu and Expa are finite, the enumeration scheme is well-defined. 

How to Find a Scheme 

First a warning: there is no guarantee that a set of patterns V possesses a finite scheme. In fact, 
based on empirical evidence, most don't (and if they did, their depth would be so large that the 
'polynomial' in the "polynomial growth guarantee" would be of such a high degree as to make it 
almost as bad as exponential growth. 

But that's what so nice about non-trivial human research. If we know beforehand that we are 
guaranteed to succeed, than it is not research, but doing chores. So our "algorithm" is not known 
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to halt. To make it a genuine algorithm, we make the maximal depth part of the input, and restrict 
the search to Prefix Schemes of bounded depth. If we fail, then the program returns 0. 



The algorithm for looking for a Prefix Scheme for a given set of patterns P is a formalization of 
the procedure described above. The input is a set of patterns, and an integer, MaximalDepth. 

We start with the empty permutation (p as belonging to Expa. Its only refinement is 1. Unless 
7^ = 0, 1 would also belong to Expa. Its refinements are 12 and 21. We examine each prefix 
permutation a, in turn, and see whether it has a non-empty J := C[a] (forced relations), and 
equipped with that J, we look for a reversely deleteable place. This we do by looking at all the 
conceivable events that the place under consideration can participate in, and look at all the implied 
events, hoping to see amongst them an event that does not include the examined place. If indeed 
each and every possible calamity in which the examined place participates in, implies another event 
in which it does not participate in, (by using the transitivity of the order relation), then that place 
is indeed reversely deleteable. 

If such a reversely deleteable place exists, then a becomes a member of Redu, and the lucky place, 
that made it a member is recorded as A[a]. If there are no such places, we reluctantly put a in 
Expa, and find its refinements, that we store as B[a]. We then examine these in turn, until we 
either get a prefix-permutation of length MaximalDepth -|- 1, in which case we sadly exit with 0, 
or all the offsprings of the member of Expa are in Expa U Redu. 

Important Remark: All the above deductions are made completely automatically by the com- 
puter. So we have (a very primitive) 'Artificial Intelligence'. 

Using the Maple package WILF 

The procedure in the package WILF, that looks for Schemes for a set of patterns, is Scheme. The 
syntax is: "Scheme(Set_Of_Patterns,MaximaLDepth)" . For example, to get a scheme for V = 
{123,132}, type: "Scheme({ [1,2,3] , [1,3,2] }, 2) ;" (without the quotes, of course), followed, 
of course, by Carriage Return. Having found the scheme (let's call it sch), to find the number of 
permutations on n objects avoiding the patterns V, you do ^Miklos{n, sch); < CR >". For example 
"Miklos (3 , sch) ; <CR>" should yield 4. To get the first L terms of the sequence a{n;V) (after 
sch := Scheme{V; Depth) was successful for some Depth), do "SchemeSequence (sch,L) ; <CR>" . 
The package WILF could also try to guess (empirically, but perhaps rigorizably) , a linear recurrence 
with polynomial coefficients. This procedure is called SchemeRecurrence. The reader should look 
up the on-line help. 

SchemeF does what Scheme does, but more generally, by trying Scheme all on the images of V 
under the dihedral group consisting of inverse, and reverse. 



This paper's website littp://www. math. temple. edu/~ zeilberg/WILF.html has numerous sam- 



ple input and output files. 
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The Maple package HERB 

This is the non-rigorous counterpart of WILF. There a scheme is found empiricaUy by checking 
(Julian) for small n, and extrapolating. It was written before the rigorous WILF, and helped a lot 
in developing the latter. 

Future Directions 

Prefix Schemes are equivalent to Suffix Schemes, but it should be possible to have mixed schemes. 
Also one may be even more creative in partitioning A^. Hopefully, a more general notion of scheme 
would be more successful. Also, it should be possible to em,pirically guess generating functions (or 
perhaps, redundant generating functions (in the sense of MacMahon)), for the A^, in ReduL) Expa, 
which, once guessed, are rigorously provable. This, in particular, would entail a constant-term 
expression for A{n; V), which, by using the WZ method, should lead to a rigorous derivation of the 
recurrence guessed by procedure SchemeRecurrence in WILF. 
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